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We study the bulk deformation properties of the Skyrme nuclear energy density functionals. 
Following simple arguments based on the leptodermous expansion and liquid drop model, we apply 
the nuclear density functional theory to assess the role of the surface symmetry energy in nuclei. To 
this end, we validate the commonly used functional parametrizations against the data on excitation 
energies of superdeformed band-heads in Hg and Pb isotopes, and fission isomers in actinide nuclei. 
After subtracting shell effects, the results of our self-consistent calculations are consistent with 
macroscopic arguments and indicate that experimental data on strongly deformed configurations 
in neutron-rich nuclei are essential for optimizing future nuclear energy density functionals. The 
resulting survey provides a useful benchmark for further theoretical improvements. Unlike in nuclei 
close to the stability valley, whose macroscopic deformability hangs on the balance of surface and 
Coulomb terms, the deformability of neutron-rich nuclei strongly depends on the surface-symmetry 
energy; hence, its proper determination is crucial for the stability of deformed phases of the neutron- 
rich matter and description of fission rates for r-process nucleosynthesis. 

PACS numbers: 21.10.Dr, 21.60.Jz, 21.65.Ef, 24.75. +i 
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I. INTRODUCTION 



Ongoing efforts to develop the nuclear Energy Den- 
sity Functional (EDF) of spectroscopic quality are faced 
with the challenge to find optimal experimental and the- 
oretical constraints that would help us pinning down its 
various coupling constants. Traditional parameteriza- 
tions of nuclear EDFs such as Skyrme and Gogny of- 
ten rely, through the fitting protocol applied, on a com- 
bination of carefully selected experimental data as well 
as pseudo data characterizing properties of nuclear mat- 
ter (see, e.g., discussion in Refs. [J-Hl)- Considering 
the fairly simple forms (and small number of coupling 
constants) of commonly used nuclear EDFs, the success 
of the nuclear density functional theory (DFT) to de- 
scribe a wide range of nuclear properties has been truly 
remarkable. However, the robustness of these parame- 
terizations when going away from the neighborhood of 
the stability valley, where experimental data are abun- 
dant, to the neutron-rich region where data are scarce or 
non-existent, is questionable. Indeed, large differences in 
predictions for very neutron-rich or super-heavy elements 
seen for various EDFs is highly unsatisfactory. In 

fact, recent systematic studies of Skyrme EDFs showed 
that some coupling constants cannot be properly con- 
strained by existing data, and that the current forms of 
EDFs are too limiting [l^l- Moreover, early attempts 
to employ statistical methods of linear-regression and er- 
ror analysis [llj have been revived recently and applied 
to determine the uncertainties of EDF parameters, er- 
rors of calculated observables, and to assess the quality 
of theoretical extrapolations [1, [l^, fl^ - fl^ . The major 



uncertainty in nuclear EDF lies in the isovector channels 
that are poorly constrained by experiment; hence, new 
data on unstable nuclei with large neutron excess hav- 
ing a lar ge "lever arm" from the valley of stability are 
essential [15l |. 

The pool of fit-observables entering the optimization 
protocol of EDF, usually contains experimental data 
characterizing both bulk (global) and local nuclear prop- 
erties, as well as theoretical pseudo-data pertaining to 
global nuclear matter properties (NMP). The charac- 
terization in terms of "bulk" and "local" is not very 
precise and somehow arbitrary; it has its origin in the 
macroscopic-microscopic approach, which offers a de- 
scription in terms of a macroscopic liquid drop (whose 
properties change smoothly as a function of nucleon num- 
bers) and shell correction that oscillates rapidly with shell 
filling [T6l - [l9j . In the context of DFT, the binding energy 
of a nucleus of mass A and neutron-excess / = [N — Z) j A 
can be split into a smooth function of / and A, and a 
fluctuating shell correction term by means of the Struti- 
nsky energy theorem [l^, [20l - [23j . This theorem, together 
with the shell-correction method, offers a formal frame- 
work to link the self-consistent DFT with macroscopic- 
microscopic models which often provide useful insights 
in terms of the liquid drop (or droplet) model and shell 
effects. 

Single-particle shell structure can be accessed through, 
e.g., experimental separation energies and single-particle 
strength. Such fit-observables are often used in the de- 
termination of EDF parameterizations P, [13] . Similarly, 
the volume and symmetry terms of the liquid drop model 
(LDM) [2^, [2^ are related to nuclear matter proper- 
ties, which are relatively well determined. This effec- 
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tivcly constrains specific combinations of parameters of 
the EDF. However, surface terms of the LDM are harder 
to pin down. In a recent work the leptodermous 
expansion of the nuclear binding energy was revisited, 
and the LDM parameters of several EDFs were extracted 
from DFT calculations for very large nuclei. While the 
surface and curvature terms came out to be fairly robust, 
it was found that the surface-symmetry coefficient agsym 
of the LDM could vary by a factor 3 across the set of 
parameterizations. This coefficient enters the expression 
for the symmetry energy in the LDM: 



LDM 
sym 
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^sym 



^ssym 
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(1) 



where agym is the (volume) symmetry energy coefficient. 
A thorough compilation of symmetry energy coefficients 
obtained in different EDFs, assuming various definitions, 
can be found in Sec. 4.7 of Ref. [28| . 

In early Hartree-Fock (HF) and extended-Thomas- 
Fermi studies [29l - [3^ using various EDFs, a correlation 
between agsym a-nd pointed out. Namely, EDFs 

having large values of asym have also large |assym|- Since 
the surface-symmetry coefficient is negative, these two 
terms act in opposite directions in i?sym- The correlation 
between bulk and surface symmetry energy was further 
discussed in Refs. f28l. Issl - fsq : it was concluded that the 
presence of the correlation makes an absolute determina- 
tion of flsym and flssym from nuclear masses difficult (see 
also discussion in an early Ref. [37|). 

The experimental information about agsym is fairly lim- 
ited. The ratio of the surface-symmetry to symmetry (or 

/asym, has 



-''ssym / Ltsym; 



volume-symmetry) coefficients, rs/v 
been estimated from the electric dipole strength distri- 
bution [ll] , masses [H, |39| - |4]| , masses and radii [sp . Issj , 
and excitation energies of isobaric analog states [42| . Re- 
cently, an attempt has been made to extract assym [43j 
from the separation energies through the displacement 
of neutron and proton chemical potentials. They noted 
a large A- variation oi rs/v- As discussed later, the DFT 
values of rg/y obtained in Refs. [13, El] are fairly 
consistent with phcnomenological estimates. 

Since the absolute value of agsym is not well con- 
strained by experimental data on ground-state (g.s.) nu- 
clear properties, one needs to find some mechanism that 
would enhance the surface-symmetry term with respect 
to the dominant volume symmetry energy. Since the 
surface-symmetry energy increases with both neutron ex- 
cess and the nuclear surface area, it is anticipated that 
strongly deformed configurations of nuclei with appre- 
ciable neutron excess can be of help. Indeed, the nuclear 
shape deformation increases the surface area thus ampli- 
fying the surface-symmetry energy in a neutron-rich nu- 
cleus. Conversely, the precise determination of surface- 
symmetry energy is important to describe the deforma- 
bility of neutron-rich systems and validate theoretical ex- 
trapolations. In this context, one can mention several 
phenomena involving neutron excess and deformation: 

• Position of the neutron drip line. Deformed nuclei 



are expected in several regions near the neutron 
drip line [4^,|4^. In some cases, deformation energy 
can impact their mere existence. For instance, it 
has been predicted that there exist particle-bound 
even-even nuclei that have at the same time nega- 
tive two-neutron separation energies due to shape 
coexistence effects |45| . 

• Borders of the superheavy region. The super- and 
hyperheavy nuclei with Z >126 can exist in states 
associated with very exotic topologies of nuclear 
density as the competition between Coulomb, sur- 
face, symmetry, and shell effects can give rise to for- 
mation of voids [13 . The subject of exotic (bubble, 
toroidal, band-like) configurations in nuclei with 
very large atomic numbers has been addressed by 
several studies (48l - [5]| . It is difficult to say at 
present whether these exotic topologies can occur 
as metastable states [H^l and what is their stability 
to various shape deformations. 

• Fusion and fission of neutron-rich nuclei. Synthesis 
of heavy and superheavy neutron-rich nuclei is pro- 
foundly affected by nuclear deformability through 
the energetics of fusion and fission valleys (52l . [53| . 
Our ability to describe fission of neutron-rich sys- 
tems is also important for modeling nuclear reac- 
tors. While DFT calculations are currently able 
to predict barrier heights of known nuclei with a 
typical accuracy of 20%, the resulting uncertain- 
ties in fission cross sections are still large |46[ . The 
ability of modern nuclear EDFs to predict neutron 
induced-fission rates for neutron rich nuclei that 
cannot be measured is crucial. At this point, the 
dependence of rates on fission barriers is apprecia- 
ble m. 

• Rotational properties of neutron rich nuclei. Nu- 
clear deformation determines the response of the 
nucleus to angular momentum. Little is known 
about the collective rotation of very neutron-rich 
systems [ssj and the corresponding interplay be- 
tween deformation, isospin, and rotation. 

• Astrophysical r-process. Fission of neutron rich nu- 
clei impacts the formation of heavy elements at the 
final stages of the r-process through the recycling 
mechanism [s^ [s^l- The fission recycling is be- 
lieved to be of particular importance during neu- 
tron star mergers where free neutrons of high den- 
sity are available [H^, [13- Also, neutrino-induced 
fission of neutron rich nuclei could affect the r- 
process flow in some scenarios [ssj . 

• Structure of neutron stars. Nucleonic phases in the 
inner crust of neutron stars are associated with very 
neutron-rich deformed nuclei or strongly deformed 
pasta and anti-pasta phases (sol . [60| . 

The primary motivation of this paper is to assess the 
role played by the surface-symmetry energy in neutron 
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rich nuclei. To this end, we study the surface-symmetry 
contribution to the LDM energy and demonstrate that it 
may be as important as the Coulomb term in driving de- 
formation properties of very neutron-rich systems. Using 
the self-consistent DFT, we survey excitation energies of 
superdeformed (SD) states in the Hg-Pb region and SD 
fission isomers in the actinidcs. By subtracting micro- 
scopic shell corrections, we extract the macroscopic part 
of the HF deformation energy and demonstrate that the 
outcome is consistent with qualitative macroscopic esti- 
mates. Our results indicate that experimental data on 
strongly deformed configurations in neutron-rich nuclei 
are key for optimizing the isospin channel of the nuclear 
energy density functional. 

Our paper is organized as follows. In Sec. |lll based 
on the spherical and deformed LDM, we discuss gen- 
eral properties of the surface-symmetry term. In par- 
ticular, we analyze symmetry-energy parameters of var- 
ious Skyrme EDFs, study global correlations between 
Ossym and asym , and show how to disentangle the surface- 
symmetry term by studying deformed configurations in 
nuclei with nonzero neutron excess. Self-consistent DFT 
calculations of fission isomers in the actinides and band- 
heads of super-deformed rotational bands in the ^ ^ 190 
region are presented in Sec. IIIIl The methodology used 
to extract the smooth contribution to the total energy is 
outlined in Sec. IIVI Section |V] presents our calculations 
of the smooth contributions to the deformation energy of 
SD states and compares HFB and LDM results. Finally, 
the conclusions of our work are given in Sec. IVII 



A. Approaches to Bulk Nuclear Properties 

Some LDM parameters arc fundamental NMPs and 
can be determined microscopically from ab-initio calcula- 
tions of the equation of state of nucleonic matter (66l - [68| . 
Another, phenomenological strategy is to obtain LDM 
constants, or at least some of them, from a direct fit to 
selected experimental data from finite nuclei. The origi- 
nal work of Myers and Swiatecki followed such a strategy 
[69j : by modeling local fluctuations in particle numbers 
due to shell effects, one can extract smooth LDM trends 
from experimental nuclear masses. Subsequent refine- 
ments involved the upgrade from a simple drop to a more 
accurate droplet model [26j], which allowed to pin down 
additional terms in the leptodermous expansion. Further 
refinements can be found in, e.g., Refs. [sol [g^. I70l - l72j . 

Just as in the microscopic approach, whose outcome 
depends on both the input (i.e., nucleon-nucleon in- 
teractions) and the theoretical method used to solve 
the many-body problem, the results of phenomenolog- 
ical procedure depend on the choice of fit-observables 
and prescription used to compute shell corrections (see 
Ref. [7^ for a recent concise overview of this topic). 
There are significant correlations among the different 
LDM terms and some parameters are poorly determined 
j36j . In particular, precise extraction of higher-order 
isospin-dependent terms requires abundant data for very 
neutron-rich and/or heavy nuclei, which are not available 
at present. 



II. LIQUID DROP MODEL FROM THE 
SKYRME ENERGY DENSITY FUNCTIONAL 

The LDM provides a convenient parametrization of the 
bulk part of the binding energy of a spherical even-even 
nucleus with Z protons and N neutrons. Expressed in 
terms of volume (flvoi), surface (csurf), curvature (ccurv), 
symmetry, surface-symmetry. Coulomb, and Coulomb 
exchange parameters, it reads: 



l2/3 



^1/3 
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5^ch^i/3 a,\2-k) 5rghyli/3' 



where e is the electric charge and rg the Wigner-Seitz 
radius. The justification of ^ can be given in terms of 
the leptodermous expansion valid for systems with a well 
developed surface |6ll - l63| that sorts the various contribu- 
tions to the binding energy of finite nuclei in terms that 
have transparent physical meaning. The expansion ^ 
can be extended to higher orders |64 [65j . The second- 
order symmetry energy term oc is not always included 
in the macroscopic LDM but it naturally present in the 
microscopic LDM expression (27j . 



B. Spherical Liquid Drop Based on Density 
Functional Theory 

An advantage of the macroscopic approach to bulk nu- 
clear properties is that it can also be applied in the con- 
text of the nuclear DFT. While some LDM constants 
pertaining to infinite or semi-infinite NMP can be ex- 
tracted from EDF parameterizations [ll, [2^ [t!] , surface 
and curvature terms are best determined by using the 
semi-classical approach [75l[76j or by removing the contri- 
bution from shell effects from self-consistent DFT results 
[23} . There are relatively few examples of latter studies in 
the literature, and most were confined to spherical sym- 
metry. In Ref. [T^I, the parameters a vol, Osym, Osurf, Ossym 
and r'^ were estimated from spherical HFB calculations 
using the finite-range Gogny force DIS [zl]. In [t^, a 
similar work was carried out for the NL3 parametriza- 
tion of the Relativistic Mean Field (RMF) jsot . 

T his p rogram was carried out more systematically in 
Ref. [27[ for Skyrme EDFs and several parameterizations 
of RMF Lagrangians using the HF method. The main 
difference with [73, [z3 ^ti^'t the convergence of the lep- 
todermous expansion was tested using an extended sam- 
ple of very large spherical nuclei. The Coulomb terms 
were ignored to be able to approach nuclei of arbitrary 
sizes and to avoid radial instabilities (Coulomb frustra- 
tion) characteristic of systems with many protons. The 
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shell corrections were extracted from the HF results ac- 
cording to Green's function method and the generalized 
Strutinsky smoothing procedure of Refs. [8ll - l83| . Table H] 
displays the values of the symmetry, surface-symmetry 
and surface coefficients for various realizations of the nu- 
clear EDF. These values are fairly close to those obtained 
in Rcf. [25. We also note that in the case of NL3, there 
arc relatively large differences between the values in Ta- 
ble H] and those reported in Ref. : these can be at- 
tributed to different ways of extracting the shell correc- 
tion, and fitting the LD formula (sample size, treatment 
of the Coulomb term). 



TABLE I: Surface, symmetry, and surface-symmetry LDM 
coefficients (in MeV) of various EDFs extracted from lep- 
todermous expansion in Ref. [27| . The Skyrme EDF 
parametrizations are: SkM* [H, SkP [H, BSkl [H, BSk6 
[131, SLy4-SLy6 SkI3-SkI4 [111 SkO The RMF 

Lagrangians are: NLl [9l[, NL-Z [ol and NL-Z2 [H. For 
comparison, the results of LDM fits are given: LDM'^' [s^ 
and LDM'2)-LSD [zil. 



EDF 




Q-sym 


flssym 


EDF 


QiSUT 


Clsym 


Qissyni 


SkM* 


17.6 


30.04 


-52 


NLl 


18.8 


43.48 


-110 


SkP 


18.2 


30.01 


-45 


NL3 


18.6 


37.40 


-86 


BSkl 


17.5 


27.81 


-36 


NL-Z 


17.8 


41.74 


-125 


BSk6 


17.3 


28.00 


-33 


NL-Z2 


17.4 


39.03 


-90 


SLy4 


18.4 


32.01 


-54 










SLy6 


17.7 


31.96 


-51 


ldm(^) 


21.1 


30.56 


-48.6 


SkI3 


18.0 


34.84 


-75 


ldm(2) 


19.4 


29.28 


-38.4 


SkM 


17.7 


29.51 


-34 


LSD 


17.0 


28.82 


-38.9 


SkO 


17.3 


31.98 


-58 











As noted in Ref. 127|, the leading surface and sym- 
metry terms appear relatively similar within each family 
of EDFs, with a clear difference for agym between non- 
relativistic and relativistic approaches. Obviously, even 
small variations in Ogurf and agym seen in TablcUcan have 
an appreciable impact on the binding energy, as these 
coefficients are multiplied by large A- and /^-dependent 
factors For the surface-symmetry coeSicient, how- 

ever, there are much larger discrepancies. For Skyrme 
EDFs, for instance, there is a factor of two between the 
largest and smallest value. This demonstrates that assym 
is very poorly constrained in the current EDF parame- 
terizations (see also discussion in (28j). 

In addition, as discussed in the previous section, there 
appears a clear (anti-)correlation between the (positive) 
value of the symmetry coefficient and the (negative) value 
of the surface-symmetry coefficient |28h36| . Figure [TJ 
displays the pairs (assym, isym) for the Skyrme EDFs of 
Tabic Hand the EDFs of Rcf. [H]. The ratio rs/v ex- 
tracted from experimental masses is rs/v ~ —1-7 [36j . 
When lookirig into details, however, it exhibits a large 
A-variation [43|; rg/y ranges between — 1 (for A > 12) 
and -1.7 (for A > 50). As discussed in [13], the data 
on the electric dipole polarizability yields rs/v * — f .65. 
The DFT values shown in Fig.[T]arc not inconsistent with 
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FIG. 1: (color-online) Correlation between the symmetry and 
surface-symmetry coeflicients extracted from Skyrme EDFs 
from Table HI (dots) and Skyrme EDFs of Ref. (circles). 
The phenomenological LDM values of Table [llare also indi- 
cated (stars) as well as the hydrodynamical [3a] and mass 



estimates: r 



s/v 



-2 and 



[43,14 

lines). The linear fit to the values of Ref. 
dotted line. 



1.6, respectively (dashed 
28| is shown by a 



these phenomenological estimates. While a correlation 
between agym and Ossym is clear, a very large spread of 
values is indicative of the inability of current data on g.s. 
nuclear properties to adequately constrain flgsym- R is in- 
teresting to note that the LDM values of Table U and phe- 



nomenological estimates cluster around Ogy 



30 MeV 



and Os 



-45McV. 



To get more insights into the consequences of this ob- 
servation, we plot in Fig. [5] the symmetry and surface- 
symmetry contributions to the binding energy per nu- 
clcon E^^^/A for the microscopic LDM derived from 
Skyrme EDFs listed in Table |T| along the LDM val- 
ley of stability. The latter one is defined by minimiz- 
ing the LDM energy in the {Z, N) plane. The symme- 
try energy slightly increases with mass number (upper 
panel), owing to the fact that the valley of stability bends 
down for increasing Z, so that the energy per nucleon 
Esym = flsym-f^ iucrcascs. Contrariwise, the surface- 
symmetry energy slightly decreases with A because of 
the A^^/"^ dependence. 

Interestingly, the total symmetry energy (i.e., the 
sum of volume and surface contributions) exhibits much 
smaller spread between various EDFs: from about OA/ A 
McV for both E'symM and Essym/A to about 0.1 /A MeV 
for the sum. This is a consequence of the aforementioned 
correlation between volume and surface symmetry ener- 
gies that implies that large discrepancies between indi- 
vidual contributions to the bulk energy tend to cancel 
out at the level of the total binding energy. We note that 
nuclear binding energies are indeed prime fit-obscrvables 
constraining parameters of most EDFs. 



> 



0.5 




1.5 



1.0 



0.5 







1 Surface symmetry energy 



— SkP ....... SLy6 — Skl4 

Skl3 — SkM* - o... SkO 
-»-BSk6 ....... SLy4 -— BSk1 



Total symmetry energy 




50 100 150 200 250 
Mass number yA 

FIG. 2: (color-online) Contributions to the microscopic LDM 
energy per nucleon along the LDM valley of stability: vol- 
ume symmetry term asym/^ (top), surface-symmetry term 
assyml^ A~^^^ (middle), and the total symmetry energy (bot- 
tom), for the microscopic LDM derived from Skyrme EDFs 
of Table m 



as nnea^ 



pnLDM TpLDM TpLDM 



{bs ~ l)a,u,fA2/3 + (6curv - l)acurvAl/3-f 



(6, - l)assym/'A2/3 



(3) 



where the geometrical factors bs, 6curv and be depend on 
the shape of the deformed drop (by definition, they are 
equal to unity at the spherical shape). Since the nu- 
clear volume is conserved in the LDM, the surface and 
curvature 6-factors increase with deformation. On the 
other hand, 6c becomes less than one as the Coulomb 
energy of the deformed drop is lower than that of the 
spherical drop. While the volume-symmetry energy is 
shape-independent, the surface-symmetry term has the 
same dependence on the nuclear shape and A as the sur- 
face term. Consequently these two contributions to the 
symmetry energy behave differently in deformed nuclear 
drops. 



100c 




C. Deformation Energy of Nuclear Liquid Drop 



The relative LDM contributions of the volume- 
symmetry and surface-symmetry energies can be disen- 
tangled if shape deformation is present. Indeed, the de- 
formation energy of the deformed LDM can be written 



FIG. 3: (color online) Individual contributions to the total 
deformation energy per nucleon of the microscopic LDM for 
nine Skyrme EDFs for two A = 100 isobars: i°°Sn (/ = 0, 
top) and ^™Zn (/ — 0.4, bottom). The assumed quadrupole 
deformation is /?2 ~ 0.6. 

As an example, we plot in Fig. [3] the individual con- 
tributions to E^^f^ for the two A = 100 drops at a 
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fixed quadrupolc deformation P2 = 0.6. (All remain- 
ing deformations are set to zero.) Specifically, shown 
arc contributions from the surface, curvature, Coulomb, 
and surface-symmetry terms corresponding to different 
Skyrme EDFs. The Coulomb radius ro was assumed to 
be the same as the Wigner-Seitz radius defining the sat- 
uration density. 

Since for -'^'^"Sn the isospin excess is zero, the defor- 
mation energy contribution coming from the surface- 
symmetry term vanishes. The variations of the total 
LDM energy between different EDF paramctrizations are 
relatively small and primarily related to slightly different 
values of agurf- The picture changes dramatically when 
going to ^™Zn, a very neutron rich nucleus with isospin 
/ = 0.4. The variations between predictions of different 
EDFs have a much larger amplitude and are caused al- 
most exclusively by the surface-symmetry term. This in- 
dicates that the theoretical differences in the LDM defor- 
mation energy of heavy neutron- rich nuclei are almost en- 
tirely driven by the poorly determined surface-symmetry 
term. As it is clear from Fig. [3J an experimental access 
to this term can be provided by extracting shell energy 
from the measured masses of very deformed configura- 
tions in neutron-rich nuclei. Another interesting obser- 
vation drawn from the deformed LDM exercise is that, 
contrary to the usual scenario in which the macroscopic 
deformability is solely driven by the competition between 
surface and Coulomb terms, the macroscopic deformation 
energy of very neutron rich nuclei involves a three-way 
competition between the repulsive surface term and at- 
tractive Coulomb and surface-symmetry terms. 



III. SKYRME HFB CALCULATIONS AT 
LARGE DEFORMATIONS 

The binding energy of a deformed nuclear configura- 
tion can be decomposed into a macroscopic part and 
shell correction. In order to determine whether the 
macroscopic features related to the surface-symmetry en- 
ergy, identified in Sec. Ill CI show up in self-consistent 
DFT calculations for actual nuclei, we have performed 
Hartree-Fock-Bogoliubov (HFB) calculations for a num- 
ber of states at large deformation with available experi- 
mental information. 



Survey of Superdeformed Bandheads and 
Fission Isomers 



We selected two regions of the nuclear chart: in the 
actinides, there is a number of isotopes where the excita- 
tion energy of the fission isomer is relatively well-known 
[osj . In the neutron-deficient Hg and Pb isotopes, the 
linking transition between the SD and g.s. bands have 
been identified for several nuclei, so that the energy of 
the 0+ band- head could be extracted [96l-[l0l|. All SD 
band-head data used in this work are listed in Table. |lTl 



TABLE IL Experimental energies of 0^ band-heads of SD 
states in A=190 mass region and in the actinides. 



Nucleus 


Esd(0+) (MeV) 


Reference 


^"^Hg 


5.3 (9) 


[98] 




6.017 


m 


192p^j 


4.011 


[96J 


194p^j 


4.643 


[99] 


196pb 


5.630(5) 


[m 


236 u 


2.750 


[95J 


238 y 


2.557 


[95J 


240p^ 


2.800 


[95J 


Cm 


1.900 


[95J 



HFB calculations were performed with the DFT solvers 
HFBTHO [102] and HFODD glOSrrlOTj . To benchmark 
binding energies of superdeformed configurations, we em- 
ployed several Skyrme EDFs in the particle-hole channel. 
Pairing correlations were modeled by an mixed-pairing 
interactio n with a dependence on the isoscalar density 
jlOSl Il09l | . All calculations were performed with a cut-off 
energy of Ecut = 60 MeV to truncate the quasi-particle 
space. For each parametrization of the Skyrme interac- 
tion, the pairing strength was fitted to the average neu- 
tron pairing gap in ^^°Sn according to the procedure de- 
scribed in Ref. |l08| . In both solvers, the quasi-particle 
solutions to the HFB problem are expanded on the de- 
formed Harmonic Oscillator (HO) basis. Since we are 
probing very elongated systems, we performed the cal- 
culations using a stretched basis with a large number of 
deformed HO shells, A'max = 20. All calculations were 
performed assuming axial, reflection symmetric shapes. 
The constrained minimization was p erformed using the 
augment ed Lagra ngian method [110| and the procedure 

of Rcfs. [ml . [112} . 

As an illustration of typical deformation landscapes 
in the two regions. Fig. U shows the calculated poten- 
tial energy curves for ^^''Pb and ^^ey functions of 
the quadrupole deformation /3 extracted from the mass 
quadrupolc moment (Q20) smd the total rms radius; 



5 (r2) ■ 



(4) 



While the actinide nuclei of interest are always predicted 
to have prolate-deformed ground states with ~ 0.3, 
neutron-deficient Hg and Pb isotopes show a more com- 
plex g.s. patte rn in volving coexisting oblate and spheri- 
cal structures [ll3l |. 

The predicted excitation energy of the SD minimum 
relative to the g.s., E^^^ ~ E^b — Eg,s., between the 
HFB SD minimum and g.s. minimum can be com- 
pared with the experimental value E* . The residuals 
AE = E;^ - E*^p are plotted in Fig. M for 22 different 
Skyrme EDFs. It is rather striking to notice that, for a 
given nucleus, the differences between various EDFs can 
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FIG. 4: (color-online) Potential energy curves for ^^''Pb (top) 
and ^^^U (bottom) versus quadrupole deformation calcu- 
lated with SkI3, SkI4, SkM*, SkO, SLy4 and SLy6 Skyrme 
EDFs. All curves are normalized to the spherical point. Ax- 
ial symmetry is assumed. 



be as high as 4 MeV, which is often greater than the 
excitation energy itself. These large fluctuations some- 
times occur within a family of Skyrme EDFs, e.g., SLy[x], 
and has been explained in some cases by the different 
recipes to treat the center of mass |ll4| . By contrast, the 
Brussels-Montreal paramctrizations Bsk[x] and Msk[x] 
are more consistent with one another. An appreciable 
EDF-dependence for SD states had already been pointed 
out in previous Refs. |ll5l . llld |. Similarly, the sensitivity 
of fissio n barriers on EDF paramctrizations was studied 
in Refs. |ll7l - lll9| . In the context of this work, it is espe- 
cially interesting to poin t ou t that the surface-symmetry 
term has been claimed (ll8l | to have a significant influ- 
ence on self-consistent fission barriers. 

The large spread in calculated values of E* comes as 
little surprise: very few EDFs have been optimized by 
considering data probing large deformations. The im- 
portance of considering strongly deformed shapes when 
fitting the co uplin g constants of EDFs was discussed in 
Refs. [H, 113, Il24l | but this program has been carried out 
only in a handful of cases. The SkM* functional [U has 
been fitted by considering the experimental fis sion bar- 
rier in Pu. The DIS Gogny interaction [zl, was 
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FIG. 5: (color-online) Residuals AE — El^^ — E*^p (top) and 
rnis deviations from experiment (bottom) for various Skyrme 
EDFs. Additional refe rences for th e Sk yrme f orces : SLy5- 
SLy7 [Mj^SkMP KM, SkX-SkXC [HJ, SIII [Hil, MSkl- 
MSk6 [l23ll. BSk2 [sj SkO' HO]. The average rms deviation 
for the set of EDFs considered (marked by a dashed line in the 
lower panel) is 1.26 MeV for the nine data points of Table [III 
and 1.34 MeV for the fission isomers nuclei. 



also optimized for fission p roper ties. In the Bskl4 EDF 
of the HFB-14 mass model [126| . data on fission barriers 
were utilized to optimize the EDF parameters by adding 
a phcnomcnological collective correction accounting for 
the zero-point rotational-vibrational motion. In this ar- 
ticle, we do not employ zero-point corrections as we are 
primarily interested in the deformation p rope rties of the 
functionals themselves. We refer, e.g., to jl27j for a more 
thorough discussion of dynamical correlations and their 
impact on deformation properties of nuclei. We note in 
passing that such correlations are supposed to impact 
standard DFT predictions of g.s. energies of Hg and Pb 
nuclei due to coexistence effects |l28l . [1291] . 
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B. Estimation of Theoretical Errors 

Since the values of AE in Fig. [5] are subject to nu- 
merical and experimental uncertainties, it is important 
to estimate their respective errors before assessing the 
model dependence of results. In order to validate E^-^, 
we studied the convergence of our HFB results with re- 
spect to the size of HO space used. Figure [HI shows the 
HFB+SkM* energy of the g.s. and SD state of ^^Opu cal- 
culated with the HFBTHO solver as a function of A^max, 
The HFBTHO numbers are compared to the benchmark 
results obtaine d wit h the precise coordinate-space DFT 
solver HFBAX poj . 

With the large HO basis used here (A^niax=20), the 
theoretical error on the energy of either the g.s. or the 
SD state is around 600 keV. Since the HFB theory is 
variational, the error on the excitation energy is in fact 
much smaller (see the bottom panel of Fig.[S|), and comes 
principally from the different convergence rates of HFB 
states with (32 ~ 0.3 and /32 ~ 0.6. Those differences arc 
due to the combination of effects coming from the basis 
deformation and the choice of oscillator frequency hu. 
At A^max ^ 16, the dependence on hu is rather weak; 
hence, the only remaining source of fluctuations is the 
basis deformation. For the residuals, we estimate the 
latter empirically to be at most 500 keV for 16 shells, 
and less than 100 keV for A^max > 16 shells. 
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FIG. 6: (color-online) Convergence of the HO basis expan- 
sion for the HFB+SkM* binding energy of fission isomer and 
ground state (top panel) and excitation energy of fission iso- 
mer (bottom panel) in ^""'Pu as a function of the HO basis 
size. Results are compared with the benchmark numbe rs ob - 
tained with the precise coordinate-space solver HFBAX [1301 ] . 



The convergence pattern of HFB calculations seen in 
Fig. |6] is to a large extent exponential. A similar be- 
havior h as also b een observed in ab-initio calculations 
of Refs. [llMll. However, in all these many-body ap- 
proaches, the size of the actual model space grows combi- 
natorially with the number of active particles and single- 
particle states taken, which is not the case for DFT. 
In a recent work, the exponential convergence of wave- 
functions expanded in a HO basis has i n fact been related 
to its weak differentiability properties |l34l |. It has been 
argued therein that this may be a generic property of 
systems with exponentially decaying wave-functions. 

For all the nuclei considered in Fig. [5] and Table [ill 
experimental g.s. masses are known to an excellent pre- 
cision of approximately 2 keV |l35| . In the A 190 
region, the uncertainty of the SD band-head comes from 
the extrapolation of the rotational band down to spin 0+ . 
This procedure is slightly mo del-d epen dent, but its error 
is estimated to be only ^^5 keV 96l - ll0l| . In the actinides, 
the determination of the excitation energy of the fission 
isomer is slightly less precise: it is about 5-10 keV for 
236,238u and about 200 keV for "^^^Pu and Cm 

Considering the above, the theoretical fluctuations of 
'^several MeV in AE seen in Fig. |S]are well above numer- 
ical uncertainties in E^^^ and experimental uncertainties 
in E*^p. Consequently, these deviations are rooted in ac- 
tual EDF parametrizations. In the next section, we shall 
investigate the relation between the fluctuations in AE 
and the underlying LDM description. 



IV. BULK DEFORMATION ENERGY OF THE 
SKYRME ENERGY DENSITY FUNCTIONAL 

In order to extract the smooth LDM energy from HFB 
results, the fluctuating contributions to the energy (i.e., 
shell-correction and pairing terms) must be removed. Af- 
ter describing the details of the extraction technique em- 
ployed, we show how the leptodermous expansion of the 
smoothed HFB energy works. 



A. Pairing and Shell Corrections 

To extract shell and pairing corrections from the total 
HFB energy is not an easy task as the building blocks of 
HFB are quasi-particles, rather than the single-particle 
states that enter the Strutinsky energy theorem [2^ . 
Moreover, while the contribution of pairing correlations 
to the total energy must be eliminated, the induced shape 
polarization must be kept, as it is relevant for making the 
direct comparison with experiment. 

To extract the effect of HFB pairing, we carried out HF 
calculations at the equilibrium deformations of HFB g.s. 
and SD configurations. This was achieved by constrain- 
ing the expectation values of the HF multipolc moments 
at those corresponding to HFB solutions. In practice, it 
sufficed to consider Q207 Qio-, Qeo, smd Qsq moments. 
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higher-order terms being neghgible. The advantages of 
this procedure are twofold. Firstly, it enables us to re- 
move all direct pairing effects. Secondly, it provides a set 
of single-particle HF states that can be used to compute 
the shell correction SEshoii- 

Shell corrections were calculated according to the pro- 
cedure outlined in Rcfs. [H, H^. It combines the stan- 
dard shell-correction smoothiiig method (our original im- 
plementation is based on Ref. [131) with the Green's func- 
tion oscillator expansion method technique that is aimed 
at removing the spurious contribution to SEghcii stem- 
ming from the non-resonant continuum of positive en- 
ergy states. Following Ref. [8l|, we employed the follow- 
ing smoothing parameters: smoothing widths 7„ = 1.66 
for neutrons and jp = 1.54 for protons (in units of 
huJo = 4:1/ A^^^) and the curvature correction p = 10. 
This choice guarantees that the generalized plateau con- 
dition is satisfied 1821. 
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FIG. 7: (color-online) Proton smoothed density g{e) in the 
g.s. and SD configuration of ^''"Pu calculated with the SkM* 
EDF and N^^^ = 16 HO shells. 



As an illustration. Fig. [7] shows the smooth density 
g{e) for protons in the g.s. and fission isomer of ^'^°Pu 
obtained with SkM* EDF. The generalized plateau con- 
dition requires that this function must be linear across 
several oscillator shells, and this is indeed well fulfilled. 
Figure [8] displays the convergence of the shell correction 
contribution to the deformation energy, iJ-EfiJ^n — i^-E'fhoii' 
as a function of A^max- While the convergence is not 
perfect, the uncertainty remains fairly small, around 200 
keV. 

We should emphasize that there exist alternative ways 
to extract shell correction, see e.g., Refs. [73l . ll36l |. How- 
ever, since we use LDM parameters extracted in Ref. [13] 
by employing the Green's function prescription (82j . it is 
important to follow the same procedure in order for our 
analysis to remain consistent. 
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FIG. 8: (color-online) Convergence of the shell correction 
contribution to the deformation energy, — 5£'fhcii5 fo'' 

^''''Pu as a function of A^'max. Calculations were performed 
with with the SkM* EDF. 



B. Determination of Microscopic LDM 
Deformations 

To compare the LDM deformation energy with the 
HFB bulk energy, one needs to properly define the shape 
of the (sharp) surface of the drop. For the axial, re- 
flection symmetric shapes considered in this work, the 
drop surface is typically parametrized in terms of de- 
formations /3; defining the multipole expansion of the 
drop radius. The problem consists therefore in map- 
ping a set {/3;}z=2,4.- - .Afq to a set of multipole moments 
(Qa)^^ = {r^Y\o) coming from HF calculations. 

The LDM deformation parameters /?/ can be d eter- 
mined from the system of non-linear equations jl37l | 



{Qxim^iQxr' A-0,2,4,...,iV, 



(5) 



However, in such an approach involving standard (vol- 
ume) multipole moments, the role of higher-order mul- 
tipoles becomes artifici ally exag gerated. It was there- 
fore argued (see Refs. [l38l . Il39l | and references quoted 
therein) that a mapping between the two sets of shape 
deformations can be best achieved by using the surface 
multipole moments defined as Qx = r^Y\Q, which have 
a much softer radial dependence than volume moments. 
Deformation parameters /3/ can therefore be extracted by 
requiring that the set of equations for dimcnsionless sur- 
face moments Il38l. Il39i 



(Qa(/30> _ (Qa 



,HF 



,2\HF 



(6) 



be satisfied, with l.X = 2, . . . , Nq. In practice only the 
three lowest terms with Z = 2, 4, 6, 8 are important at SD 
shapes. This choice provides the best mapping between 
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first extract the spherical charge radius Rfj 



^1/3 



from the self-consistent spherical HF calculations for each 
nucleus {Z,N) and then use this value of Tq'^ in Eq. ([3]). 
Since the charge radius thus obtained does not come 
from a systematic fit but is obtained locally, it intro- 
duces shell fluctuations into LDM. However, since spher- 
ical self-consis tent radii behave smoothly as a function 
of shell filling jl40| the corresponding shell effect is very 
small indeed. 

Our actual determination of the LDM charge radius 
goes as follows. From the spherical rms proton HF radius 
{Rp), we extract the rms charge radius according 
to the standard formula: 



(7) 



where (r^) = 0.764 fm^ and (r^J = — 0.116fm^ are the 
squared rms proton and neutron charge radius, respec- 
tively. The geometrical charge radius Rq^ is then ob- 
tained from the rms charge radius in the usual way: 



ch 



(8) 



We note that the condition (i?^h)^° = (-^ch)"^ together 
with Eq. ^ implies that the charge surface of the mi- 
croscopic liquid drop is close to that of HF. 



D. Example: Extraction of the Smooth 
Deformation energy for 236u 

To illustrate the extraction procedure of smooth defor- 
mation energy from HF results, and to assess the quality 
of the leptodcrmous expansion, we carry out constrained 
HF+SLy4 calculations for ^^^U shown in Fig. O The 
constraint on the quadrupolc moment was determined 



> 



self-consistent multipole moments and deformations of 
the sharp surface. 



C. Coulomb Polarization 

In Ref. [131, the microscopic LDM parameters were 
extracted from a set of spherical HF calculations with- 
out the Coulomb term. However, such a methodology 
is clearly not applicable in realistic calculations. First, 
the Coulomb term crucially affects nuclear deformabil- 
ity. Second, while its contribution to the total energy 
can easily be subtracted, the Coulomb potential induces 
a long-range polarization of the mean-field, which af- 
fects the equilibrium deformations, single-particle states, 
etc. Most importantly, this Coulomb polarization is 
deformation-dependent. As a result, the contribution of 



the Coulomb term to the excitation energy, -EJouj >2 



vary by be up to several MeV for the interactions that 
we consider in this study. 

To take this effect into account at the LDM level, we 
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FIG. 9: (color-online) Extraction of the LDM deformation 
energy from constrained HF+SLy4 calculations for ^-^^U at 
several values of quadrupole deformation /3. Shown are: the 
total shell correction (5i5shcii (squares), smooth HFB defor- 
mation energy E^ciil^) (dots), and the corresponding LDM 
deformation energy 

EdS^W) (triangles). The inset shows 
the equivalent LDM deformations Pi with Z=4,6, and 8. 



so that the deformation (3 of Eq.(j4]) takes the values 
/3 = -0.20, -0.15, . . . , +0.80. Since the HF potential en- 
ergy curve consists of several sharply-crossing configura- 
tions as evidenced by rapidly varying LDM deformations 
shown in the inset of Fig. [HI we made no attempt to inter- 
polate between the mesh points in /3. The shell correction 
was extracted at each f3 according to the procedure dis- 
cussed in Sec. IIV Al The smooth energy at deformation 
/3 is given by: E^^{/3) = E^^{l3) - SEshdiW), which 
defines the smooth component of the HF deformation 
energy Ef^f{P) = E^^{P) - E^^{P = 0). At each p, 
the LDM deformation parameters /3; are computed ac- 
cording to Sec. IIVBI Finally, the deformed LDM energy 
E^^f^{/3) is obtained using Eq. ^ with LDM constants 
taken from Table U and the charge radius defined accord- 
ing to Sec. IIV CI As expected, the smooth deformation 
energy is growing as a function of deformation; the local 
variations arc due to configuration changes in HF calcu- 
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lations and corresponding changes in higher-order shape 
deformations. 

Even though the agreement between E^^^ and E^^f 
is not perfect, it is gratifying to sec that the LDM energy 
nicely follows the smooth energy extracted from HFB. 
The deviation has multiple sources such as the error on 
the regression analysis carried in [27| , uncertainties of the 
shell-correction procedure, neglect of the second-order ef- 
fects in density fluctuations [22| , LDM assumption of the 
sharp surface, limitations of the leptodermous expansion 
used, etc. Considering all this, the quality of the lep- 
todermous expansion of the deformation energy in de- 
formed nuclei is very reasonable. 



V. SURFACE-SYMMETRY ENERGY AND 
DEFORMED NEUTRON-RICH NUCLEI IN DFT 

We are now ready to determine the smooth part of 
the deformation energy from HFB results and compare 
it with the microscopic LDM using the methodology de- 
scribed in Sec. IIVI In Fig. |3] we demonstrated that the 
uncertainties in the determination of Cssym impact the 
deformation properties of neutron-rich nuclear drops. It 
is therefore interesting to see whether in the studies of 
well deformed states in realistic nuclei these uncertain- 
ties would show up. 
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FIG. 10: (color-online) Smooth excitation energy E* of fis- 
sion isomers in 236,238^^ 240p^^ g^^^ 242 calculated in HFB 
and LDM for seven EDFs. See text for details. 

FigurefTUlshows the smooth excitation energy E* of fis- 
sion isomers calculated in HFB and LDM for seven EDFs. 
For each nucleus, we first carried out HFB calculations 
to determine the g.s. and SD configurations. The con- 
strained HF calculations are then performed based on the 
multipole moments of the HFB solution. Shell energies 
are subtracted from g.s. and SD HF energies, and this 
defines the smooth part of the excitation energy E^-p^ in 
HFB, i.e., the smooth deformation energy of the excited 



state relative to the ground-state. Using the surface mo- 
ments obtained in the g.s. and SD minima of HFB, we 
extract the equivalent LDM deformation parameters /3/ 
and the LDM excitation energy. 
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FIG. 11: (color-online) Same as in Fig. [in] except for SD 
bandheads in W2,i94Hg and wa.w^Fb. 

As in Sec. IIIIl clear differences between various EDF 
paramcterizations can be seen. Overall, these variations 
can be as large as 4 MeV at the LDM level. As dis- 
cussed earlier in Sec. IIVDI there is a ^2 MeV shift of 
the LDM curves with respect to the HFB results. How- 
ever, it is rewarding to see that the shift is systematic 
and the EDF-variations seen in HFB are properly cap- 
tured by the equivalent LDM. The results for the SD 
bandheads are displayed in Fig. [TT] In these lighter nu- 
clei, the agreement between HFB and equivalent LDM is 
better on the average, but local fiuctuations can be ap- 
preciable (see SkO or SkI3 results for Pb isotopes) and 
might be related to a complex pattern of g.s. equilibrium 
deformations in these nuclei. 

The results shown in Figs. [10] and [11] combined with 
the overall picture of the residuals in Fig. [SJ demonstrate 
that large differences between Skyrmc EDFs exist when 
it comes to deformation properties of nuclei. While these 
differences certainly depend on variations of EDF param- 
eters controlling the shell structure, such as, e.g. the ef- 
fective mass or spin-orbit splitting, our analysis indicates 
that there are also fundamental discrepancies at the level 
of the bulk energy. One may therefore question whether 
EDF optimization protocols based exclusively on a small 
amount of data in nuclear matter and spherical nuclei are 
able to capture the deformability of EDF. 

As discussed in Sec. [IT] and in particular in Figs. [T] 
and [21 ftssym varies very significantly from one EDF to 
another. Consequently, surface and symmetry properties 
of EDFs are intertwined in a non-trivial way when it 
comes to deformability. Guided by the results of Figs. [TUj 
and 1111 we may wonder whether the large variations in 
flsswri are indeed reflected in the results of self-consistent 
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FIG. 12: (color-online) Surface and surface-symmetry con- 
tributions to the LDM excitation energy of fission isomers in 
the actinides compared to the smooth HFB excitation energy 
for the same Skyrme EDFs as in Figs. 1101 and [TT] All curves 
are shown relative to SkM* results. 
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FIG. 13: Relative contributions of the Coulomb, surface- 
symmetry, curvature and surface terms to the equivalent LDM 
excitation energy of SD states in a sequence of U isotopes. 
Calculations are based on SkM* and BSk6. 



calculations. 

Figure [T^ shows the surface and surface symmetry con- 
tributions to the LDM excitation energy of SD states in 
the actinides for the same Skyrme EDFs as in Figs. [10] 
and 111! The equilibrium deformations that are used in 
the LDM for both the ground-state and SD state are ob- 
tained in HFB. The LDM results are compared to the 
smooth HFB energy E^-p^. To facilitate interpretation, 
all curves are normalized to SkM* values. In this way, we 
can better compare relative variations obtained in vari- 
ous EDFs. It is interesting to see that the inter-EDF 
fluctuations of i^HFB rather well correlated with the 
surface-symmetry energy. In other words, the contribu- 
tion from the Coulomb and curvature terms (not plotted 
in Fig. [T2] for better legibility) cancel out the surface 
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FIG. 14: Quadrupole deformations /3 in SkM* for ground 
states and fission isomers for the same U isotopes as in Fig. 
1131 (top) and the corresponding excitation energies (bottom) . 



term to a large extent. This result is significant because 
it seems to confirm the simple analysis of Sec. [IT] in a re- 
alistic case: in nuclei having large neutron excess / (here 
of the order of / « 0.2), differences in deformation energy 
between various EDF paramctcrizations reflect the differ- 
ences of the surface-symmetry coefficient. By contrast, a 
similar analysis of individual macroscopic contributions 
in the Hg-Pb region does not allow to pin down a single 
LDM term as a primary deformation driver. 

To further illustrate the importance of the surface- 
symmetry term, we calculated the LDM excitation en- 
ergy of the fission isomer for a sequence of U isotopes. 
Here, employed the SkM* and BSk6 parametrizations. 
The SkM* E DF i s known to perform rather well for fis- 
sion barriers jl4ll Il42| . Its surface-symmetry coefficient 
is also close to the average among Skyrme forces and val- 
ues from phenomenological estimates, so it can be viewed 
as fairly representative of the Skyrme functionals. The 
BSk6 parameterization gives a reasonable rms deviation 
for excitation energies of SD states, see Fig. [5] The iso- 
topes considered include some very neutron-rich species 
important in the context of r-process fission recycling. 
It is worth noting that HFB potential energy landscapes 
change considerably within this isotopic sequence. For 
example, A^=184 is a neutron magic number for SkM* 
[Slf : hence, g.s. configurations around ^^^U are spherical 
rather than prolate, see Fig. [T31 The equivalent LDM 
equilibrium deformations, therefore, reflect these struc- 
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tural changes in a non-trivial way. 

For each isotope, we computed the relative contribu- 
tion of the surface, surface-symmetry, curvature, and 
Coulomb term to the total LDM excitation energy: these 
are the only terms that depend on deformation. Fig- 
ure [13] shows the percentage of the LDM excitation en- 
ergy from these contributions. As the neutron excess 
grows, one can notice the gradual relative decrease of the 
Coulomb contribution - which depends only on the num- 
ber of protons and proton density, and therefore remains 
relatively constant in value. This decrease is compen- 
sated by an increase of the surface-symmetry contribu- 
tion. For the most neutron-rich nuclei considered here, 
the total contribution from i?ssym is as large as 19% for 
SkM*. For the two paramctrizations selected here, the 
role of the surface-symmetry term increases by a factor 
3 from 236 u ^t^^, fission recycling region. 



VI. CONCLUSIONS 

This work contains a comprehensive study of defor- 
mation properties of nuclear energy density function- 
als based on the Icptodcrmous expansion of the smooth 
nuclear energy. Since symmetry and surface-symmetry 
terms in the expansion arc strongly correlated, a way to 
resolve them is to consider data on deformed neutron- 
rich nuclei, in which the surface-symmetry term is am- 
plified. Based on intuitive LDM arguments, we argue 
that deformation properties of neutron-rich nuclear drops 
are governed by an interplay of the deformation-driving 
Coulomb and surface-symmetry terms, and the surface 
energy that acts against shape deformation. To estimate 
this interplay, we extracted the smooth deformation part 
of the HFB energy by means of the shell correction pro- 
cedure. 

Self-consistent DFT calculations for excitation ener- 
gies of SD states show marked differences in their predic- 
tions depending on the parametrization used. For the set 
of EDFs considered, the average rms deviation between 
predicted energies of SD states and experimental values 
is 1.26 MeV. Within this set, the MSkl parametrization 
provides the best overall reproduction of the data: the 
corresponding rms deviation is 0.53 MeV, and this set a 
benchmark for future improvements. For the subset of 
fission isomer data, the best performer is SkM: its rms 
deviation is 0.48 MeV. 

We showed that inter-parametrization differences re- 
flect to a large extent macroscopic properties of EDFs. 
In particular, our calculations indicate that the bulk 
deformation properties of actinides are strongly driven 
by surface-symmetry effects, while in the proton-rich 
A ~ 190 nuclei there is more competition between the 
various macroscopic contributions. This finding should 
have an impact on the fissility of heavy, very neutron- 
rich nuclei of the kind encountered e.g. in the r-process. 
For example, the surface-symmetry contribution to the 
bulk part of the excitation energy of fission isomer in 



very neutron-rich uranium isotopes can reach ~20% as 
compared with '-^5% for ^■^^U. 
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FIG. 15: The upper part of the chart of the nuclides with the 
X = I limit indicated for SkI 3 and SkM* EDFs, the value of 
77 = 1.7826 used in Ref. [l43l ] (LDM), and no isospin depen- 
dence (»7 = 0). The region of known nuclides is marked by 
black squares. 

The importance of the surface-symmetry term on fis- 
sion can be quantified at the LDM level by the dimen- 
sionless fissility parameter: 



-Ecoui(sph) 

2£;,urf(sph) ^ 47^1(1-77/2)' 



(9) 



where 77 



i/isurf- If a:: > 1, the nuclear liquid 



drop is unstable to fission. In the presence of neutron 
excess, the fissility parameter increases, i. e., th e tendency 
towards fission increases. In Refs. [5^, Il43l |. the value 
77 = 1.7826 was used. By taking LDM parameters from 
TableUwe see that 77 is 1.9 for BSk6, 2.9 for SkM*, and 
4.16 for SkI3, i.e., this parameter is very uncertain. 

Figure [T51 shows the LDM fission limit for the SkI3 and 
SkM* EDFs, as weU as for 77 = 1.7826, and 77 = (no 
isospin dependence). The minimum value obtained for 
Z^/A, i.e., 47 in Eq. is not very precise as it de- 



pends on assumptions about the LDM constants jl44 |. 
Therefore this diagram should be considered as a quali- 
tative guidance. A clear message drawn from Fig. [15] is 
that the surface-symmetry term can significantly impact 
LDM fission barriers: the greater the value of 77, the lower 
the threshold for fission. This result is especially impor- 
tant in the context of the fission recycling mechanism in 
the r-process and hot fission reactions leading to excited 
neutron-rich superheavy nuclei. Since shell effects are t o 
a large extent washed out at high temperatures [ll, ll45l | , 
the fission of hot compound nuclei is expected to be gov- 
erned by the LDM fission barrier (or smooth HFB de- 
formation energy) . As seen in Fig. [T5] the uncertainty 
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in flssyin, hence rj, makes it difficult to reliably predict 
fission rates of the heaviest and superheavy neutron-rich 
nuclei. (In this conte xt we note that according to the 
recent estimates |l46l | assym and r] are expected to very 
weakly depend on temperature.) 

The results obtained in this paper suggest that adding 
to the list of fit-observables data on strongly deformed 
nuclear states (such as excitation energies of SD states 
or fission barriers), combined with the usual constraints 
on bulk properties and shell structure, should constrain 
quite effectively the surface properties of the nuclear 
EDF. Such a strate gy is currently being pursued within 
the UNEDF project jlivl 1 148j . On the experimental side, 
new information on deformed properties on neutron-rich 
systems is the key. 
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